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Abstract 

The aim of this short review is to give an introduction to inonopoles and to present theoretical 
derivation of two particular monopoles in ferromagnetic metals, a hedgehog monopole and a spin 
damping monopole. In an electromagnetism in the vacuum, described by Maxwell's equations, 
magnetic field and electric field are not symmetric, since there is no monopole, a particle having 
a finite magnetic charge. Still, monopole has been an exciting object for a long time, and was 
discussed on a phenomenological ground by Dirac in 1931. A theoretical possibility of monopole 
generation was first given by 't Hooft and Polyakov in 1974 in a context of a symmetry breaking in 
a grand unified theory (GUT). A GUT monopole has not been discovered in experiments so far. In 
contrast to in the vacuum, several kinds of monopoles are expected to emerge in solids associated 
with various symmetry breaking mechanisms. Of particular interest is metallic ferromagnetic 
systems, since a breaking of a symmetry of conduction electron spin, described by an SU(2) algebra, 
can give rise to monopoles. Indeed two monopoles are theoretically predicted in ferromagnets; one 
is a hedgehog monopole arising from a topological spin structure and the other arising from spin 
damping in the presence of spin-orbit interaction. In this paper, we focus on these monopoles, 
while other objects similar to monopoles but not coupled to effective electromagnetic fields, such 
as spin ice monopoles, are touched only briefly in the introduction. Those monopoles are extended 
objects coupled to effective electromagnetic fields, which are described by Maxwell's equations 
with monopole contribution. The effective fields are the one coupled to spin of a particle like 
electron, muon and neutron; two monopoles are thus detectable by electric measurements. Spin 
damping monopoles can be generated in simple systems such as a junction of a ferromagnet and a 
heavy element with strong spin-orbit interaction such as Pt. This monopole is essential in coupling 
electronics with magnetism, and is thus expected to play an essential role in spintronics. 
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FIG. 1. Left: Schematic figure of a magnet and magnetic field around it. Right: A monopole, 
particle having only N (or S) magnetic pole, and the magnetic field it emits. 

I. INTRODUCTION 

A. Monopole 

Macroscopic magnets have two poles called N (north) and S (south) at two different edges 
(Fig. [1]). We also know that we cannot extract only N or S pole from magnets; when we 
divide a magnet, a pair of N and S poles is created at the edge and each piece becomes a 
magnet with the equal amount of N and S poles. There is therefore no monopole at least 
in the energy scale we concern in our life. This is because the magnets are made of spin, 
a quantum magnet the electron carries, which creates the divergencelss magnetic field only. 
In other words, each spin contain N and S poles, which are not separable since elementary 
particles are pointlike. 



Possibility of a particle carrying only N or S pole, called magnetic monopole, was first 
discussed by Dirac [1|. He showed based on a macroscopic consideration that monopole is 
tied to a string of singularity, and its existence is allowed only if the charge of the monopole is 
quantized by unit of 2Tih/e (e is the electron charge) . A microscopic mechanism of monopole 
creation was theoretically discovered by 't Hooft [2| and Polyakov |3| independently. They 
showed that monopole arises in a grand unified theory (GUT) of elementally particles when 
the symmetry breaking of the SU(5) symmetry to the U(l) symmetry of the electromagnetic 
field occurs. Such monopoles may have been created in the early universe (at about 10~^° 
sec after the big bang) and may still be around us. So far, however, no evidence has been 
obtained in experiments waiting for a monopole from universe to go through superconducting 



detectors 



B. 



J] or detecting the ionization [5|. The energy scale of the GUT monopole is 10 



17 



GeV, and so creating one in accelerator on the earth is impossible. 

B. Maxwell's equations without monopole 

Before starting to discuss monopoles, let us discuss the law of conventional electromag- 
netism without monopole. We know that electric field, E, and magnetic field, B, behave 
differently, i.e., electric field has finite divergence while magnetic field is divengenceless |6|. 
To put in equations, 

W.E=P 
e 

W-B = 0, (1) 

where p is charge density and e is dielectric constant. The two fields are governed by another 
set of equations describing rotational components; 
^ ^ dB 

W xB = pj + e/i— , (2) 

where /x is magnetic permeability and j is current density. The fields E and B are therefore 
not symmetric; magnetic field is generated by electric current by the Ampere's law, while 
there is no driving current for electric field. One may imagine that the electromagnetism law 
is more beautiful if there is a current that drives electric field via an analog of the Ampere's 
law. Let us write such current as jm- The first equation of Eq. ([2]) then becomes 

VxE = -j„, - — . (3) 

Taking the divergence of this equation, we have 

= -V-j^-|(V-S), (4) 

which indicates that the second equation of Eq. ([T]) needs to be modified as 

V ■ B = p^, (5) 

where pm is a quantity satisfying a conservation law of 

Pm + V ■ jm = 0. (6) 





FIG. 2. Left) The Ampere's law allows us to generate magnetic field from the electric current, j. 
In contrast, there is no Ampere's law for the electric field in the electromagnetism if U(l) symmetry 
is exact. The exception is the case with monopole: monopole current, jm, generates E (Right). 

Equation ([5]) indicates that pm emits a magnetic field like a N or S pole; hence pm is a density 
of monopole. The current j^ is then a monopole current. Thus, if monopole exists, E and 
B become symmetric, and the electric field can be generated by applying the monopole 
current. 



The aim of this short review is to give an introductory description of monopoles in 
condensed matter and explain in detail monopoles in ferromagnetic metals. The paper is 
organized as follows. Monopole is discussed from a macroscopic viewpoint in ^JTTl Section IIIII 
is a brief introduction of monopoles in solids. Other topological objects similar to monopoles 
are introduced in §IV[ These objects are touched only briefly, since we would like in this paper 
to focus on monopoles which really couple to electromagnetism via Maxwell's equations. 
Sections |V] and |Vl] are the main content of the paper presenting theoretical description 
of hedgehog monopole and spin damping monopole, respectively. Relativistic notation is 
summarized in ^M and spin damping is briefly described in 



II. MONOPOLE FROM MACROSCOPIC VIEWPOINT : DIRAC'S STRING 

Let us go on to discuss monopoles from a macroscopic viewpoint, i.e., based on the 
Maxwell's equations. As is well-known, in the electromagnetism without monopole, the 
magnetic field B is expressed as a rotation of a vector potential, namely, as S = V x A. 
When monopole is present, this is no longer true, but still one can define a vector potential 



as a line integral over a certain line L starting from monopole to the infinity. Let us consider 
a case of a single monopole at r = 0, 

B{r) = ^-^, (7) 

where g is a. monopole charge. We can define a "vector potential" 

^^\r) = Y.e,,J dr'^B,{r-r'\ (8) 

where Bk is the kth component of the magnetic field (including monopole) and L is a string 
connecting r = and the infinity. The "magnetic field" of this vector potential A^^^ then 
reads 

V X A^^\r) = B{r) + j dr' p„,{r - r'), (9) 

where pm is the monopole density, Pm{r) = gS^{r) if a single monopole case. The first term 
of Eq. ([9]) correctly reproduces the field B but we have an additional non local singular 
field represented by the second term. The singular magnetic field represented by the second 
term of Eq. ([9]) exists along the string L. Dirac argued therefore that monopole needs to 
have such string (called the Dirac's string) but this string must not been "seen" by physical 
particles [l|. If many monopoles exist, many strings L each attached to a single monopole 
are needed to describe the vector potential. 

The unobservability of the string is ensured if the monopole charge g is quantized to 
certain values. One way to derive the quantization condition is to impose the condition that 



Is 



the shift of the Dirac's string does not modify the wave function of the charged particles [6|. 
Equation ([8]), explicitely written in a single monopole case reads 

^'(r) - iE^.4 *^^- (10) 

jk -^^ I I 

If one shift the string L to be another string L', the vector potential changes to 

jk -"^ I I 

where C = L' — L is a closed path surrounded by L and L' . The last term of the right-hand 
side is proportional to the derivative of the solid angle f2c(^) subtended by C observed at r. 




FIG. 3. When a solid angle of a closed path C observed at P is 0,c, it is 47r — f^c* at P' in the 
opposite side of the path. 

In fact, a difference of solid angles observed at r and r + 6r, which we call 6Qc', is (noting 
that dr' x 6r is a vector normal to the plane having a length of an element of area) 

5nc= [ (dr' X 6r) ■ -^—^. (12) 



c 



We therefore see that 



\r — r 



(r — r' 



VQc = dr'x } ^, (13) 



c 



\r — r 



and that Eq. flTT|) becomes 
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^iL') ^ ^(L) ^ JLvi^c;. (14) 

Ait 

Since a phase for the electron's wave function is given by an integral of the vector potential, 
this modification of the vector potential results in a phase change of e~^^'^. Quantum 
mechanics requires that this phase is single-valued. However, solid angle Qc is not single- 
valued. Multi-valuedness of Qc occurs when Qc changes by Air as an electron goes through 
the plane spanned by C (Fig. [3]). The single- valuedness of the electron's phase is thus 
ensured by requiring that 

eg = 27rnh, (15) 

where n is an integer. This is the quantization condition for the monopole charge, pointed 
out by Dirac. 

The above argument is based on the assumption that the magnetic field is written as a 
rotation of a U(l) vector potential even in the presence of monopoles. This assumption, 

6 



however, is wrong, since the exact U(l) gauge invariance does not exist when monopole 
exists. As we will show later, the monopole field in the case of the symmetry breaking 
mechanism is expressed by the additional field orthogonal to the U(l) degrees of freedom. 
The Dirac's string containing a singular magnetic field (the last term of Eq. (jH])), is thus an 
artifact arising when one tries to describe the monopole field within the U(l) gauge theory 
assuming B = V x A, and it is eliminated when one takes account of microscopic origin of 
the monopole. 

In the case of 't Hooft-Polyakov monopole, monopole is created from the symmetry 
breaking [7| . Their argument is for a system of non-Abelian gauge field coupled to a Higgs 
represented by a vector field, 0° {a = x, y, z). They showed that there is a solution where 
Higgs field behaves at the infinity as 

0° oc — (r ^ oo). (16) 

r 

In this configuration, magnetic field is shown to contain monopole, and monopole density is 
determined by Higgs field as 

Pm fx ^ ^ eijueap-^di(t)'^dj(t)^dk(P. (17) 

ijk aj3^ 

A volume integral of this monopole density is written as a surface integral at the infinity, 
and it turns out to take only integer value called a winding number. This topological nature 
of the Higgs field leads to a quantization of monopole charge. Thus, 't Hooft-Polyakov 
monopole emerges from a symmetry breaking from non-Abelian group to U(l) group driven 
by a condensation of Higgs field. Dirac's string is not necessary here, since monopole field 
is created from Higgs field, which is different from a gauge field. 

In ^|Vl we will discuss in detail a monopole arising from a symmetry breaking in ferro- 
magnets. 

III. MONOPOLES IN SOLIDS 

A. Hedgehog monopole 

As we have mentioned, in the electromagnetism in the vacuum, monopole predicted at 
high energy has not been found. Even if it is found in the future, we cannot make a device 






FIG. 4. Schematic illustration of a magnetization structure for a hedgehog monopole, discussed 
by Volovik. This structure has a singularity at the center, and is topologically non-trivial. The 
electron coupled to this magnetization structure feels the effective electromagnetic field with a 
monopole. 

since its energy is too high. In contrast to electromagnetism in the vacuum, we have hope 
in solids. In fact, electrons in solids feel another U(l) gauge field which couples to the 



B 



electron's spin |8|. The spins are object in an SU(2) space, and thus if the breaking of its 
symmetry occurs, a resulting effective U(l) gauge field may contain monopoles according to 
the 't Hooft and Polyakov scenario. In addition, the energy scale of the symmetry breaking 
is at low energy, less than leV, and so device application would be straightforward. 

Such a monopole was indeed theoretically pointed out to exist in ferromagnetic metals 
by Volovik in 1987 |9|. In fact, a monopole arises from a strong sd coupling between the 
conduction electron and local spin, which specifies the projection of the conduction electron 
spin having SU(2) symmetry to a U(l) space. When local spin structure is a topologi- 
cally non-trivial one called a hedgehog, the projected effective electromagnetic field contains 
monopoles. 

B. Spin damping monopole 



Recently, another monopole in ferromagnetic metals was predicted 10|, lUj . The idea is 
to include spin-orbit interaction besides sd interaction. Spin-orbit interaction modifies the 
projection to the U(l) plane defined by the sd interaction, and thus new monopole may arise 

8 




FIG. 5. Schematic illustration of a spin damping monopole. A damping of spin (S) generates an 
orbital angular momentum L for electrons when spin-orbit interaction (SO) acts. A circular motion 
due to the orbital angular momentum is equivalent to that caused by an electric field induced by 
a monopole current jm via Ampere's law. 

when the spin-orbit interaction is included. The monopole generation in such case cannot be 

discussed by use of gauge fields, since the spin-orbit interaction is not a gauge interaction. 

nstead, novel method based on a transport calculation was applied by Takeuchi and Tatara 



10| . They derived effective electromagnetic fields by calculating the electric current by use 
of Keldysh Green's functions, and showed that the fields satisfy Maxwell's equation with 
monopole. In their unique approach, a knowledge of gauge invariance was not necessary 
to explore the structure of the electromagnetism. This fact may sound surprising, but 
is natural, since a U(l) gauge invariance is equivalent to charge conservation law, which 
is strictly observed in transport calculations. In their case, a deviation from exact U(l) 
invariance due to spin-orbit interaction resulted in a monopole contribution. Based on the 
transport method, they revealed that monopole arises when the spin structure is dynamic 
and when spin-orbit interaction is included. More specifically, monopole arises when there 
is a damping of spin, represented by a damping vector (see ^ |B|) . 

N = SxS, (18) 

(S being local spin), and thus the monopole was named the spin damping monopole. 

A physical mechanism of spin damping monopole is understood as a transfer of the spin 
angular momentum to orbital one (Fig. |5]). In fact, a vector aN, where a is the Gilbert 
damping parameter (proportional to the strength of the spin-orbit interaction if in metals), 
represents the spin angular momentum dissipated. This lost angular momentum is converted 
into the orbital motion of the electrons by the spin-orbit interaction, inducing the circular 



^ 



FIG. 6. A scalar spin chirality subtended by three spins represents non-coplanarity. Spin chirality 
reduces to spin Berry's phase in the slowly varying limit. 

orbital motion of the electrons. This circular motion is regarded as a result of a fictitious 
magnetic field due to a monopole. /^From this intuitive explanation, the monopole density 
should appear when the spin damping is spatially inhomogeneous. This is indeed consistent 
with the result of Ref. 10| , which showed that the monopole density is pm oc V- A/". Similarly, 
when the spin damping is time- dependent, a monopole current arises since a temporal change 
of effective magnetic field is equivalent to a monopole current. The spin damping monopole 
was argued to be essential in the spin-charge conversion in spintronics phenomena such as 



the inverse spin Hall effect 12 1 



These monopoles in ferromagnetic systems are discussed in detail in [|V]and fVT 



IV. RELATED TOPOLOGICAL OBJECTS IN SOLIDS 

A. Spin chirality 

There are several objects closely related to monopoles. Spin chirality is one example. In 
a spin system on a lattice, a non-coplanarity of three spins. Si, S2 and S3, is represented 
by a scalar product called a scalar chirality, x = Si ■ {S2 x S3) (Fig. |6]). When the spin 
structure is slowly varying in space, spin chirality can be expanded as 

Six) ■ iS{xi) X 5(x2)) ~ J2^xi - x)^{x2 - x),S{x) ■ (d^Six) X d^Six)), (19) 



/lU 
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FIG. 7. Schematic illustration of a magnetic vortex (left) and a skyrmion (right). These are 
two-dimensional structures of magnetization with topological meaning. The topological number of 
a vortex is the number of spin rotation as one travels along the perimeter, called a vorticity. In a 
skyrmion structure, the spins at the perimeter is fixed to be parallel (spin down in this figure), and 
its topological number is the number of spin rotation as one travels from the origin to the infinity 
along the radial direction. Vortex and skyrmion in this figure have a topological number of one. 

and the spin chirality reduces to the spin Berry phase. In the continuum limit, therefore, a 
scalar chirality is represend by a vector 

X. = 5^e,,fc5-(V,-5x Vfc5), (20) 

where the direction i is orthogonal to a plane the three spins lie. This expression of the scalar 
chirality is proportional to spin Berry's phase or a local effective magnetic field generated 
by a hedgehog monopole (see Eq. (146|) ). One should note, however, that finite local spin 
chirality does not necessarily mean that a monopole exists; monopole exisitence requires a 
surface integral of the chirality, J dS ■ x [dS is an element of surface integral), to be finite, 
and this is realized only for a three-dimensional hedgehog structure shown in Fig. HJ 

Even when monopoles do not exit, a finite spin chirality leads to interesting physics, such 
as inducing chirality-driven exotic anomalous Hall effect [13|, |lj]. When spin structure such 
as domain wall is dynamic, the effective electric field acting on the electron's spin, the spin 
motive force, is induced (see Eq. ( I46l) ). Spin motive force was detected in the case of moving 
domain wall [Iq . 

A monopole was also found in systems of conventional anomalous Hall effect driven by 
spin-orbit interaction [l6|, [l7(]. This monopole is a singularity in the momentum space and 
is not coupled to the electromagnetism by Maxwell's equations. 

11 
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FIG. 8. Left: Ground state spin configuration in a frustrated magnet with a pyrochlore lattice. 
All the tetrahedra contains two in-spins and two out-spins, satisfying ice rule. Right: An excited 
state of a spin ice. Due to a flip of a spin connecting the two tetrahedra, left tetrahedra has three 
out-spins and the right tetrahedra has only one out-spin. This configuration contains effective 
monopoles; positive and negative ones on the left and right tetrahedra, respectively. 



B. Vortex and skyrmion 



Magnetic vortex and skyrmion in magnets are spin structures having topological feature 
in two-dimensions (Fig. [7j) [18[. They have finite value of local spin chirality. A vortex is a 
structure where the number of spin rotation as one travels along the circle at the infinity (or 
the edge of the disk) is the topological number called a vorticity, which is an integer. Vortex 
state is stable in a small circular disk of submicron size since the magnetostatic energy at 
the edge is the lowest for the vortex state 19|. Switching of a vortex core by an electric 
current is a hot issue from a viewpoint of fundamental science and application to non- volatile 
memories. At present, rather large current density of 3.5 x lO^^A/m^ is necessary for a core 
flip, and the switching time is not very fast, about 20ns [20 1. 



Skyrmion was originally proposed in high energy physics. In condensed matter, a magne- 
tization structure shown in Fig. [7] is called a skyrmion 2l|. Such a structure was predicted 
to arise when the inversion symmetry is broken. Skyrmions forming a lattice were observed 



in MnSi 



22 



23[. Vortices and skyrmions have finite effective magnetic field, Bs, but these 
topological structures are nothing to do with monopole (p^ and j^ vanish), since they are 
two-dimensional objects. 
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C. Spin ice monopole 



Monopoles in frustrated magnets with pyrochlore lattice is a hot subject recently 2J]. In 
the ground state of pyrochlore spin system, each tetrahedra has two spins pointing to the 
center of tetrahedra (in-spin) and two spins pointing outwards (out-spin), satisfying so called 
an ice rule. When the system is excited, ice rule is broken and some of tetrahedra contain 
three in-spins and there arise the same number of tetrahedra with three out-spins (Fig. 
8]). In such excited states, tetrahedra with broken ice rule have nonvanishing divergence of 
magnetization, V ■ M = pM- In terms of magnetic field H = —B — M, the excited states 
of spin ice systems thus have monopoles, i.e., V ■ i? = —pM 25|]. The spin confi gura tion 
corresponding to monopoles has been observed by neutron scattering experiment 26||. In 
addition, an effective magnetic charge of a spin ice monopole was measured in Dy2Ti207 
. Spin ice system is also created recently on an artificial square lattice 29 1. 



27 



One should understand that spin ice monopole is not a real monopole, since V ■ -B = is 
strictly true even in the presence of any spin configurations; it is a monopole of an artificial 
magnetic field 2J|. Nevertheless, an spin ice monopole interpretation is highly useful to 
describe excited states in frustrated spin systems from the viewpoint of spin liquid. 



D. Monopole in topological insulator 



Recent study revealed that monopole arises at the surface of topological insulators as a 
result of a image magnetic charge when an electric charge is close to a surface according 



to the following mechanism 30 



311 • ^ surface of topological insulator is described by a 



massless Dirac Hamiltonian in (2+1) dimensions. In this system, a parity anomaly arises 
from a ultra violet divergence, resulting in a Hall current perpendicular to the applied field 



J 



I'M 



2h\m\ 



Y.e^^'d^A,, 



(21) 



^lvp 



where A represents an vector potential and ?7i is a topological mass 32|. An electric charge 
emitting a radial electric field outside a topological insulator thus creates a circulating electric 
current on the surface. A magnetic field generated by this current is equivalent to the one 
emitted by a magnetic charge inside topological insulator. Therefore, electric charge in a 
proximity with a topological insulator is coupled with a mirror magnetic charge (monopole). 
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forming a state called a dyon. This coupling of electric and magnetic charges is explained 
also by an effective Lagrangian of a topological insulator in three space dimensions (called 



a 6 term) [30 1 



Le = ^ fd'xE-B, (22) 

where 6' is a constant. This Lagrangian, derived by integrating out the electrons, indicates 
that electric field and magnetic field are coupled directly, i.e., there is a magneto-electric 
effect. Lagrangian Le is written as a surface integral, and reduces to a Lagrangian which 
describes a parity anomaly in (2+1) dimensions. 

V. HEDGEHOG MONOPOLE IN FERROMAGNETS 

A. Gauge field representation of ferromagnetic metals 

In this section, we derive a hedgehog monopole in a metallic ferromagnet in a strongly 
spin-polarized case. There are conduction electrons and local spins (magnetization). Lo- 
cal spin we consider is a classical vector field, represented by S{r,t), which depends on 
space coordinate, r, and time, t. We use a field (second-quantized) representation, where 
c = {cj^,c^y (lower index ± denotes spin and t stands for transpose) and c^ represent the 
annihilation and creation operators for a conduction electron, respectively. The free part of 
conduction electron Hamiltonian is 

i/o = /rfVct(-|-V^-e.)c, (23) 

where m is the electron mass and ep is the Fermi energy. In metallic ferromagnets, con- 
duction electrons are spin-polarized by local spin, S, due to a coupling of sd-tjpe given 
by 

H,d = -J fd^rS ■ {c^cTc), (24) 

where J is a coupling constant, and [c^crc) is the spin density of the electron (cr = ax, Cy, a^ 
represents Pauli matrix). The total Lagrangian of conduction electrons, defined as L = 
Jd^rihc'c — H, where H = Hq -\- Hsd, is 



d-^r 



( fi} 

ihc^c — I iVcp — epc^c I + Mn ■ (c^ac) 

\2m 
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(25) 



where M = JS and n = S / S is a unit vector. In ferromagnetic metals, spin polarization of 
conduction electron satisfies Mr/h ^ 1, where r is elastic lifetime of conduction electron. 



This is the adiabatic condition for disordered metals [a, ui3 



3J]. In this limit, a local gauge 



transformation to choose the electron spin quantization axis along S{r, t) at each point is 



useful 



33|. The deviation from perfect adiabaticity is then described by an SU(2) gauge 



field, which is small and we treat it perturbatively. A new electron operator a = {a+,a-Y 
is defined as 



c{r,t) = U{r,t)a{r,t), 
where [/ is a 2 x 2 matrix which we further define as 



(26) 



U{r,t)=m-(T, (27) 

m being a real three-component unit vector we will define later. The matrix satisfies U^ = 1, 
i.e., U{r,t)~^ = U{r,t). A derivative of an operator c reads 

d,c{r, t) = U{r, t) [d, + f/(r, tr'd,U{r, t))a = U{r, t) [8^ + zA,,^)a, (28) 

where a gauge field (represented by a 2 x 2 matrix) is defined as 

A,,^ = -iU{r,t)-^d^U{r,t). (29) 



In terms of spin components, A^ is written as [33 1 
^s,/. = (m X 9^m) ■ cr = ^ ^"^o^a- 

a 

By the above gauge transform, the electron spin is transformed to be 
U'^aU = 2m{m-(T) -a. 



(30) 



(31) 



The aim of our gauge transform is to let this spin to be along z-axis, i.e., U ^{n- cr)U = az 
This is satisfied if we choose 



m 



sm - cos 0, sm - sm 0, cos - , , 
2 2 2 ' 



(32) 



where [6, (p) are the polar coordinates of S. The gauge field is then obtained in a matrix 
notation with respect to spin index as 

//15 \ 






' —di^i_9 sin (p — sin 9 cos (pd^cp * 



d^6 cos — sin 6* sin 4>dfj,(, 



cos 9)d^(f) 



(33) 
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The electron part of Lagrangian is written in terms of a-electron as 



d^r 



h^ 



,t. 



iha^d |VaP + epa'a 

2m 



Ma^ 



CTM 



^2 f^2 



In the matrix notation of the spin, 



d^r 



(a^+, a^_ 



ihdt 



2m 2m. 






h^ 



where 



4± 



+ (a^+,aL) 



a; ± ^^11, 



-M + ij^Al.Vi - hAl 



s,0 



2m 



^^4".V. - /iA 



2m s,« 



s,0 






(34) 




(,35) 



(36) 



and V =V — V acts only to the field operators. 

The adiabatic limit is defined as M — )■ oo. In this limit, the minority spin electron 
has infinitely high energy and thus does not exist. Off-diagonal elements of Eq. ( l35l) are 
accordingly neglected, and the system reduces to an electron interacting with a U(l) gauge 
field, v4g^^, described by a Lagrangian ((A^)^ and (A^)^ act only as a potential) 



^ad 



d^ra). 



ihdt 



2m 



V^ + e^-M + i-— 4",V* 

/?7l 



2m 



[A^Y - hA 



s,0 



a+. (37) 



B. Hedgehog monopole arising from non-adiabaticity 

When M is finite, perpendicular fluctuation represented by A^^ exists, and there is a 
finite deviation from the U(l) symmetry. These components appears in the U(l) space as a 
singular magnetic structure, i.e., a monopole. This can be shown as follows. (The notation in 
this subsection is a relativistic one [7| , and upper and lower indices have different meanings. 
See Sec. |A]for details.) A field strength of the SU(2) gauge fields is 



a=x,y,z 

where a-component is given as 



(38) 



I3l=x,y,z 



(39) 
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where e^/s^ is the asymmetric tensor in three-dimensions. By definition, the field strength 
satisfies the following identity called the Bianchi identity: 

e^^P'^D.F.^p, = 0, (40) 

where e^vpa is the asymmetric tensor in four-dimensions and 

{D,F,,p^r ^ 9.F°^, + ^ e^p,AlFZ^,, (41) 

represents the covariant derivative of field strength. 

When away from the perfect adiabatic limit, A^^ = A^^ need to be taken account of. 
Nevertheless, when the non-adiabaticity is weak, only the z component of the field strength 
is essential, which reads 

where 

<^^"' = {A^'^A^^'y - 4^'^' 4^'^), (43) 

is an anomalous field strength representing a trace of the SU(2) gauge field. In terms of a 
unit vector n, it reads 



<^^"' = -n-{d^'nxd''n). (44) 

In the adiabatic limit, F^"'"' with spin component a = x,y are suppressed and the z com- 
ponent of the Bianchi identity ( BT]) reduces to 



e'"^^'^9,Fs,,.,, = e^^f'^d^id.Al, - d^Ai^ + $,.) = 0. (45) 

^From a gauge invariance, the effective electric and magnetic fields are defined as 

E,,i = F.fli = -Vi^o + dtAli = --n-{hx V^n) 

5s,i = (V X A'^)i = ^Y1 ^'^^^ ■ (^^'^ ^ ^'^^)- (^^) 

jk 

The yU = component of Eq. (145|) then becomes 

-V-Ss + $^e^ifc9i<|),fc = 0, (47) 

ijk 
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namely 

V ■ Bs = ph, (48) 

where 

Ph = ^ etjkdi^jk = 4 X] (^tjkldiU ■ {djTi x dku)]. (49) 

ijk ijk 

Thus the effective magnetic field has a finite divergence, i.e., a finite monopole density. 
Similarly, fi = i component of Eq. fl45p reads 

(V X E,)i + dt ■ B,^i = -jh,i, (50) 

where 

jh,i = Yl eijkidt^jk - 2dj^ok) = 7 X] (^ijkl'Ti ■ {djU X dku)]. (51) 

jk jk 

Consistency of Eqs. ( H8|) ( 150|1 is guaranteed by a conservation law for monopole, 

Ph + V ■ jh = 0. (52) 

Therefore, ferromagnetic metals having a singular hedgehog spin structure in the adiabatic 
limit contains monopole. 

C. Quantization condition of hedgehog monopole 

An important feature of hedgehog monopole is that its density and current, j^ and phi 
vanish when the length of local spin is constant. This is easily seen by noting that a unit 
vector n is described by two independent angles, 6 and 0, and that three vectors diU, djTi 
and dkn in Eqs. f H9|) and fl5T]) cannot be independent. Nevertheless, the volume integral of 
the monopole density is finite due to the surface contribution if the local spin has a hedgehog 
structure shown in Fig. |H Let us define a normalized spin gauge field including a copling 
constant ^, i.e., as (upper suffix N means north) 



A^^^£il-cose)d,<p. (53) 



The effective magnetic field of A reads 



jk 
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The total magnetic flux is then 

dS-B^ = gn, (55) 

meaning that there are n monopoles with a charge g. Here we used a relation 

J2 / (^Sie.jkTi ■ (Vj-n X Vfcn) = 2 J^ / dS^e^jk sin ^(V/)(Vfc0) = Svm, (56) 

where n is an integer. This is understood by noticing that sinOdOdcj) is an element of are 
of a sphere spanned by 6 and (p, and thus an integral "^^-f, J dSieijk sin 9(V j9){'V k(p) is a 
solid angle of a sphere (47r) times a winding number of spin structure, an integer n. This 
relation is a result of a fact that the total solid angle subtended by a spin structure is Ann, 
or equivalently a fact that the left-hand side of Eq. f l56|) is a winding number multiplied 
by Svr. The hedgehog monopole from the spin structure is thus a topological object having 
vanishing local densities. The one shown in Fig. IHhas a topological number of ra = 1. 

A quantization of monopole charge g is discussed by requiring that a gauge fleld covering 
the whole space without singularity is constructed by patching together locally deflned gauge 
flelds. In fact, Eq. fl53l) is not defined at 6' = vr (south pole), since d(f) can not be defined 
there. We can define a gauge field regular at the south pole as 

Al = -£-il + cose)d,<l>. (57) 

This field has a singularity at the north pole 6 = 0, but represents the same magnetic field 
as Eq. ( l54l) . We can define therefore only a gauge field with a singularity, if we try to 
describe the whole space by a single gauge field. The singularity is regarded as a Dirac's 
string. Instead, we can cover the whole space by patching two gauge fields, A^ and Eq. (!53|l . 
which we call A^. They are related by a gauge transformation 

A^^=Al + ^^^Q-%e, (58) 

where 

e = e-^'^'t', (59) 

is a gauge transform function. This function must be single valued, i.e., is invariant under 
— 7- + 2n. Thus a condition 

eg = 2nh, (60) 
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is imposed, which is a Dirac's quantization condition (Eq. ( !T5|) ). 

Mechanism of generation of a hedgehog monopole is essentially the same as a GUT 
monopole, namely, a symmetry breaking of a non-Abelian gauge field (SU(2) for a hedgehog 
and SU(5) for GUT). There is, however, a difference in the two models. A GUT monopole 
is a composite object of a gauge field and a Higgs field, and the monopole solution is solved 
in the same footing as a Higgs field. In a case of hedgehog monopole, in contrast, local spin 
structure is treated as a background field which is treated as not to be affected by gauge 
field. 

The effective magnetic field of a hedgehog monopole configuration, Eq. ( 146|) . is sometimes 
called the spin Berry's phase, and the electric field is known as the spin motive force |8|. 
The magnetic field has been observed for example in the anomalous Hall effect [17|, and the 
electric field has been observed by inducing magnetization dynamics such as domain wall 
motion 15|. These experiments do not, however, mean existence of hedgehog monopole. 



since as explained in gTVl 



D. Hedgehog monopole in weak exchange coupling regime 

Until very recently, a hedgehog monopole has been discussed exclusively in the adiabatic 
(strong sd coupling) regime. As we have discussed in the earlier sections, a symmetry 
breaking of spin SU(2) space to a U(l) space of electromagnetism is clearly defined in 
this regime. However, hedgehog monopole is not an object restricted to adiabatic regime. 
In fact, introduction of an exchange interaction to rotationally invariant spin results in a 
symmetry breaking even if the interaction is weak. In this section, we will demonstrate 
that the hedgehog raonopole emerges even if the exchange coupling J is small following the 



analysis in Ref. 35|. We cannot approach monopoles in this regime by a standard gauge 



field argument. Here we apply instead a novel method pointed out in Ref. [10|,|35| based on a 
transport calculation. In a transport method, effective electromagnetic fields are calculated 
by evaluating electric charge density and current density induced by local spin structure 
by use of Keldysh Green's functions. A monopole field is then identified by deriving the 
Maxwell's equation for the effective fields. 

We consider a disordered metal, and take account of the spin-independent impurity scat- 
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tering represented as 

H^= f £rc^ViC, (61) 



Vi being the impurity potential. In the following calculation, the impurities are approximated 
as random point scatterers and the averaging is carried out as 

{v,{r)v,{r')), = n,uf6%r-r'), (62) 

where n, and ui are the impurity concentration and the strength of the scattering, respec- 



tively 36|. The impurities give rise to an elastic lifetime for the electron, r, which is 



calculated as 

n 



(63) 



(z/ is the density of states per volume). The total Hamiltonian discussed in this section is 
H = Ho + H, + H,^. (64) 

The electric charge density is 

p{r,t) = -etT{c\r,t)c{r,t)), (65) 

where the bracket represents the quantum expectation value and tr is a trace over spin 
indices, and electric current density is given as 

Ji{r,t) = —{Vr -Vr')^tTG<{r,t■,r',t)l,^^, (66) 

where 

G<Ar,t;r\t')^'-{cl,{r',t')c,{r,t)), (67) 

s and s' are spin indices) is the lesser component of the non-equilibrium Green's function 



37(1 . This Green's function defined on the Keldysh contour C satisfies the Dyson's equation, 

Gss'{r, t; r', t') =6s,s'9s{r, t; r\ t') 

+ fd'r" f dt"g,ir,t-r"X) 

X [6s,s"V,{r")-JSir",t")-(Tss"] 

xGs>,Ar",t";r',t'), (68) 
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FIG. 9. Diagrammatic representation of the electric current, j, induced by sd interaction with 
local spin, S. Solid lines represent the conduction electron's Green's function and wavy lines are the 
interection with S. The first diagram is a contribution of an effective electric and magnetic fields 
local in space, and the second diagram, containing a diffusion ladder (vertex corrections) denoted 
by the gray shaded oval, results in a diffusive current (a gradient of electric charge density). 

where Gss'ir,t;r',t') = —{i/h)(Tc[csir,t)cl,{r',t')]) (Tc is the path-ordering operator on 
C) and g denotes free Green's function. This equation is solved by iteration. Here, we 
assume slowly varying magnetization profile Sg^a in space and time: the spatially smooth 
magnetization structure compared to the electron mean free path i, q£ '^ 1 (g is a wave 
number of magnetization texture), and the sufficiently slow dynamics of magnetization, 
Qr ^ 1 (r2 is a frequency of magnetization dynamics). 

To see a hedgehog monopole, it is enough to discuss the electric current to the third-order 
in the exchange coupling, J. This contribution is diagrammatically shown in Fig. |9] and is 
calculated as 



ji{r,t) 



eJ3 



irmV 



iQ-r+iQt 



30m 



q,{Q ■ q')lm{9l)' 



^q,uj ' \^q',uj' ^ '-'i 



Q—q—q',Q—uj—uj' 



) 



k,q,q' ,Q U},ui' ,fl 



^r'uqlKl' 



DV^p{r,t), 



(69) 



where 



9k 



h?k'^ 



ih 
2t 



(70) 



is the advanced Green's function [ep is the Fermi energy), Im means taking an imaginary 
component, V is system volume, and D = 2ei?r/3?n denotes a diffusion constant. The last 
term is a diffusive contribution arising from vertex corrections shown in the right diagram 
of Fig. and electric charge density, p, is calculated as 

^4 



P 



AeuJ^T 



nm V.(5.(5xV.5)^^ 



(71) 
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where (■ ■ ■)d is an average including electron diffusion, which satisfies 

(-DV' + d,){F{r,t))^ = -Fir,t), (72) 

r 

[F is an arbitrary function depending on space r and time t). 

Summing over the wave vectors in Eq. ( I69l) . we obtain the electric current, 



pfi u T Apv T t 



As was pointed our in Ref. 35|], the effective electric and magnetic fields (£?h and B\i) are 



read from the above result by comparing it with a general expression, 

3 = (l//i) V X Bh + (T,Sh - DVp, (74) 

where /x is magnetic permeability. 



2 I ^ V ( hkp 



a^ = e nr/m = — — r, (75) 

3 V \ m J 

is electric conductivity. The result is 

Eh,. = -27h5-(5x V.5), (76) 

Bt,r = l^e,,kS ■ (VjS X VkS), (77) 

where 

>.6fi!0, (78) 

and we defined magnetic permeability as 
1 1 e"^ u fh " 



Obviously, these effective fields satisfy the Faraday's law and the Gauss's law with magnetic 
monopole, 

V X ^h + -Bh = - jh 

V ■ Sh = ph, (80) 

where magnetic monopole contributions are 

Jh,i = -S^heijkS ■ (Vj5 X Vfc5), (81) 

Ph = ThCiifcV.S ■ (VjS X Vfc5). (82) 
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There is therefore hedgehog monopole in the weak sd couphng case, too, although the 
coefficients in Eq. fl82|) differ from the strong couphng hmit. It is notable that the structure 
of the electromagnetism, i.e., U(l) gauge field, (and further, that with monopole in the 
present case) is embedded in the electron transport phenomena. 

We note here that the definition of E^^, B^^, /i according to Eq. ([71]) and the dielectric 
constant e according to the Gauss's law has arbitrariness. In fact, the condition imposed 
by the transport properties {j and p) is not sufficient to fix the two effective fields uniquely, 
and an additional condition seems to be required. Nevertheless, existence of monopole holds 
true; both of the monopole current and density cannot be deleted at the same time by 



redefining the fields. Further, as noted in Ref. 35|], the product e/i is invariant; a physical 



quantity of the velocity of the effective topological electromagnetic field is uniquely given as 



h 



t'top = —^= = ■^^-^=j^—^. (83) 



If elastic mean free path, i = ^^^£J^^ is loA, epT/h = 3.3 by choosing kp^^ = 1.5A, the speed 
of the topological electromagnetic wave is thus rather large, about 900m/s. The topological 
electromagnetic field may be useful to transport spin information in a different manner from 
magnon transports. 

VI. SPIN DAMPING MONOPOLE 

As we have seen, the hedgehog monopole is a topological object and has locally vanishing 
density and current only. The hedgehog monopole hence does not locally coupled to the 
electromagnetism. Very recently, a novel monopole in magnets with locally finite density 



and current density was discovered by Takeuchi et al [10|. Such a monopole creates a 
rotational electric field via the Ampere's law (Fig. [2]), and thus it acts as an anomalous 
angular momentum source which induces rotational motion of electric charge. To realize 
such a monopole, Takeuchi et al. included the spin-orbit interaction. 

The spin-orbit interaction exists in any elements including magnetic ones, and is particu- 
larly strong in heavy elements such as platinum and gold, and at the interfaces in junctions 



where the inversion symmetry is broken 38|]. Two types of the spin-orbit interaction were 



thus considered in Re 



10|. The ffist is the one from a uniform field, E-^, namely the 



Rashba interaction 39|]. Such a field is realized at the interfaces and surfaces. The Rashba 
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interaction in metallic films recently turned out to be particularly useful for current-driven 
magnetization switching 40|,|4l|. The second is the one from a random potential, Vi, induced 
by heavy impurities. This random heavy impurity model would also simulate the effect of 
spin-orbit interaction in heavy pure metals. The total spin-orbit interaction thus reads 



Hso = -\ Jd'rc^ [(Ar^r - \,Vvi) ■ (p X (t)] c, 



(84) 



where p is the electron's momentum and A is a spin-orbit coupling constant (the subscript 
R and i characterize Rashba and impurity-induced ones, respectively). The interaction with 
the magnetization is described by Hgd- The Hamiltonian of the present system is, therefore, 
given as 



H = Ho + Hi + Hsd + Hso- 



(85) 



The electric current, j, generated in the system by magnetization dynamics is calculated by 
evaluating a quantum field theoretical expectation value of electron velocity operator. 



V = -{ih/m)V + (1//1)(ArSr - XiVvi) x a. 



The electric current thus reads 



j{r,t) =etr 



2m 



z[Ar^r - XiVviir)] X o- G<ir,t;r\t) 



(86) 



(87) 



This expectation value is evaluated by solving the following Dyson's equation, 

Gss'{r, t; r\ t') =5s,s'gs{r, t; r', t') 

dV f dt"gs{r,t;r",t") 

x{6s,s"Viir")-JSir",t")-ass" 

+ 2{[Ar^r - XiVViir")] X Vr"} ■ (Tss") 

xGs"Ar",t";r',t'). 



(88) 



This equation is solved treating A and J perturabatively to the linear and second orders, 
respectively. Contributions to current is represented by Feynman diagrams in Figs. [TD] and 

m 

We consider sufficiently slow dynamics of magnetization, namely Qt ^ 1 (fi is a frequency 
of magnetization dynamics), and assume that the magnetization structure varies smoothly 
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FIG. 10. Diagrammatic representations of electric current pumped by magnetization dynamics 
and the Rashba interaction. Solid lines represent the conducting electron Green's functions, and 
double dashed and wavy lines represent the Rashba spin-orbit interaction (E^) and the interaction 
with localized spin (S), respectively. The first three contributions correspond to effective electro- 
magnetic fields, and the last two contributions, containing diffusion ladders (vertex corrections) 
denoted by the gray shaded ovals. 



in the space compared to the electron mean free path i, i.e., g^ <^ 1 (g is a wave number of 
magnetization profile). The leading contribution in this case turns out to be 



3{r,t) 



2nV 



\~^ \~^ g-i(qi+<5f2)-r+i{ni+f^2)tj^ 



k,k' ,qi,q2 uj,Qi,Q2 



r^ART 



+7^ — ^(gi+92)£fei^fcri(?fc'ii 

-DVp{r,t). 



a|2 



(89) 



The last term is the diffusive contribution arising from the vertex corrections, where the 
electric charge density p is 



P 






V ■ (^R X (5 X 5))d. 



(90) 



Here (■ ■ ■)d represents the average including the electron diffusion, defined in Eq. (!72l) . 
Summing over the wave vectors and frequencies in Eq. (I89l) . the electric current is obtained 
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FIG. 11. Contribution to the electric current driven by a precession of localized spin and ran- 
dom impurity-induced spin-orbit interaction. Filled circles are spin-orbit interactions arising from 
random impurities (SOI), filled diamonds are nonmagnetic impurity scatterings (imp), and dotted 
lines linking filled circles to filled diamonds represent impurity average. 



as 



16ei/Ai J^e^r^ 

J = :7T7. V X [S X S) 



3ff 



Erx{Sx S)- DVp. 



(91) 



This result is rewritten by use of effective electric and magnetic fields, ^s and Bg, as 

J = -V X Bs + CTc^s - DVp, 
where the effective fields are defined as 



(92) 



^s = -aREK X N, 
B, = -PiN. 



(93) 



Here 



N = SxS, 



(94) 
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FIG. 12. Left: Schematic illustration of damping of a precessing local spin S. A vector product 
of a time derivative of spin, S, and S is a vector N = S x S. Its average is N, which is along 
the precession axis, and it represents a time-averaged dissipated spin magnitude. Right: Damping 
torque, N, results in a damping of spin precession when an external magnetic field B is applied. 



is a vector representing the spin damping torque (see ^ |Bl) (Fig. [T2ll 33l . |42| . Coefficients Or 
and (3i are 

4ez/ARj2r2 



ttR 



O-rfl^ 



A = ^, . (95) 

The effective fields calculated here are the ones acting on the electronic spin in the same 
manner as the effective fields from the hedgehog monopole. Clearly, the fields [Eq. f l93|) ] do 
not satisfy the Faraday's law and the Gauss's law of the conventional electromagnetism, but 
the ones with monopole contribution, 

V • Ss = pm, (96) 

where the monopole current and monopole density read 

jm = orV X {En xN) + AiV, (97) 

and 

Pm = -A V ■ N. (98) 

We have thus proved that a monopole emerges when spin damping occurs, namely we 
have a spin damping monopole. The spin damping monopole is a composite object made 
from a magnetization configuration in the same manner as the hedgehog monopole. The 
monopole satisfies the conservation law, Pm + V • jm = 0. 
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A. Remarks on uniqueness of effective fields 

Equation (p6|l apparently contains an arbitrariness. In fact, one may think that a trans- 
formation Eg ^ Es + (l/cTc) V X C and Bg — )■ B^ — fiC, where C is an arbitrary vector field, 
is possible without changing Eq. (1961) . However, such a transform is not allowed because of 
a gauge invariance in the original space with a higher symmetry, as is known in the case of 
hedgehog monopole, where an SU(2) gauge invariance forbids this arbitrariness. 

Validness of our definition of Eg and Bs in Eq. (!93l) is supported by the following 
argument. In Eq. fl92l) . electron diffusion (D'Vp) induces the effective electric and magnetic 
polarizations Pg and M^, respectively. By using the relation 

i-DV^ + dt)p = -a,V ■ Eg, (99) 

the electric current (eq. fl92p ) is described as the rotation of the magnetic field and the time 
derivative of the electric field, 

j = V X ifs - A, (100) 

where H^ and £>= are the fields defined as 



H, = -B, - Ms 



1 

A = eE, + Ps, (101) 

with e = —(JcT being the permittivity, respectively. In Eq. f llOip . Eg and Bg are naturally 
identified with the electric and magnetic fields, respectively. By use of these fields, two of 
the Maxwell's equations become 

V ■ i/s = Pm. (102) 

Here j^ is redefined monopole current, which is non-local. 

B. Spin damping monopole generation in a ferro-normal junction 

The spin damping monopole is unique since it does not require a particular non-coplanar 
spin structure like a hedgehog, and so it exists quite generally in magnetic systems. The 
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FIG. 13. Schematic illustration of monopole pumping and detection in a thin ferromagnetic 
film attached on a non-magnetic layer. Magnetization (S) precession is induced by applying an 
oscillating magnetic field. The amplitude of presession is represented by an angle 9. The Rashba 
field, -Er, exists at the interface and creates the monopole current, jm, near the interface. The width 
of the monopole current distribution, d, is comparable to the decay length of the magnetization 
at the interface. The monopole current induces an electric current, j, via the Ampere's law at 
the interface. The impurity spin-orbit interaction directly induces positive (-I-) and negative (— ) 
monopole charge distribution, p^, at the two edges. 

simplest candidate for creating the monopole would be a thin ferromagnetic film put on 
a non-magnetic insulator or metal shown in Fig. [131 We choose the z axis perpendicular 



to the film. T 
z direction 138 



le Rashba-type spin-orbit field would then arise at the interface along the 
We excite the precession of the uniform magnetization by applying the 



alternating magnetic field in the yz plane in the presence of static field along the x axis (fer- 



romagnetic resonance 



42|). The precession results in the spin damping vector with a finite 



time average, N, along x direction. In the present case with the uniform magnetization, 
spatial derivatives in Eqs. (197]) and ( p8|) arise at the interface and at the edges, where the 
magnetization vanishes. The Rashba interaction contributes to the DC monopole current 
at the interface as 



Jm.x 



oz a 



(103) 



where d is the spatial scale of the magnetization decay at the interface. The monopole 
current driven by the random spin-orbit impurities, on the other hand, vanishes when time- 
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averaged. The total DC monopole current thus reads 

J^=-e^^EnN, (104) 

a 

{e^ represents the unit vector along the x direction). This monopole current at the inter- 
face generates the electromotive force along the y direction via the Ampere's law for the 
monopole. The monopole density induced by the random spin-orbit interaction arises at the 
edge of the ferromagnetic film since V • A?" ^ dN^/dx is finite there. The induced monopole 
density at the two edges is 

p^=T^N, (105) 

where the sign is positive on one side of the edge and negative on the other side. The 
monopoles then produces a magnetic field along the x direction as 

K = -e,P,N. (106) 

This field creates the electric current in the y direction via the conventional Ampere's law. 
The averaged electric current density generated by the spin damping monopole [Eq. (!92l) ] 
thus reduces to 

3=-ey{a,anEn + ^)N. (107) 

We have thus found that when magnetization precession occurs in a metallic ferro-normal 
junction, a current of voltage arises perpendicular to the junction and a precession axis. As 
far as electric detection concerns, a monopole effect appears qualitatively the same as an 



inverse spin Hall effect (12|, which is a widely-used experimental method for spin current 
detection. We will examine this important point later in WIIAI 

Let us look into the monopole signal quantitatively. We define energy scales of Rashba 
and impurity spin-orbit interactions as {kp is Fermi wavelength) 

Ar = — kp vr 
e 

Ai = Xikp^Vi, (108) 

where wr = j-Eji is a potential energy due to Rashba electric field, £^r. Coefficients or 
and /3i then read 



2 



"-""*(£ (fm^- '-' 
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/3i and fxekp have dimension of Tesla-sec. Here we used u = kp'^/ep, n = kp^ /{Gt^'^) and 
(Jc = e^nr/m. 

When the spin damping arises from the magnetization precession with the frequency Q 
and the angle 9 (Fig. [13]), the magnitude of a monopole- induced current density, Eq. (11071) . 
reads 



IjI =4ekln{^] sin^ 



' ^Ar 4 1 Ai 



(110) 



Ep 3 kpd Vi 

We consider a disordered ferromagnets with J/ep ~ 0.1, e-pr/h ~ 10 and k^^ ~ 2 k{eF = 
1.51 X 10~^^J= 0.95eV and r = 7.0 x lO^^^s). Rashba interaction strength is chosen 



2m 



as Ar/£:f ~ 0.1, considering an enhancement on the surfaces and the interfaces 43[. We 
assume for simphcity that t^^ is the same order of magnitude as —^. When 6 = 30° and 

^ •' kpd VI ° ep 

VL = 1 GHz, the electric current density is thus 

IJI = 1.6 X 10^ A/m^ (111) 

which is sufficiently large for experimental detection. In addition to DC, there is AC com- 
ponent in Eq. (1921) which would be accessible by a time- resolved measurement. 
The current corresponds to an effective electric field, Eq. (jHSD, with a magnitude 

For the above values of parameters, the field is 390 V/m. The voltage induced by monopole 
when a sample width is lOOnm is thus 39/iV. The magnitude of effective magnetic field, Eq. 
(l93D,is 

\Bs\=P,n = ffiekpn(^^^ (^) , (113) 

which is estimated as 5.4m Tesla if we choose /i = 10^ x /io (/io = 47r x lO'^H/m is permeability 
in vacuum) as in the case of permeability of permalloy for the ordinary magnetic field. We 
note that /i here is a permeability for effective magnetic field and can be different from the 
one for ordinary magnetic field. For correct quantitative estimate, magnetic properties for 
the effective magnetic field needs to be investigated further. 
The monopole current, Eq. (1971) . is estimated to be 

Ij^I^^ -4x10^1 V/ml (114) 
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This value, however, cannot be converted into usual dimension of current density, A/m^, 
since permeability of the effective field is not known. 

We note that the electric current estimated here is an initial current that arises when the 
pumping of monopoles starts. When the monopole current is pumped steadily, the monopole 
accumulation grows at the edges of the system, inducing a diffusive current. The steady 
monopole distribution is then determined by the balance of this backward diffusion and the 
pumped monopole current. 

VII. DISCUSSION 

In this paper, we have succeeded in proving existence of monopoles by deriving Maxwell's 
equation by a transport calculation in §VDI and §VI[ There remains, however, a few points 
to be clarified in the transport approach to spin damping monopole. 

First is a question whether spin damping monopole has a topological meaning or not. In 
a case of hedgehog monopole, it was a topological object in three space dimensions, while 
topological nature seems to be lacking in the expression of spin damping monopole density 
(Eq. fl98|) ). What is crucially different in spin damping monopole is that it needs dynamic 
spins (S). Thus, if there is any topological meaning, it should be discussed in both space 
and time. 

Second point is a gauge field representation of spin damping monopole. In a case of 
hedgehog monopole, it was obvious in the strong coupling limit. Even in the weak coupling 
limit, it is straightforward to see from the effective fields (Eq. (177|) ) that the gauge field 
describing hedgehog monopole is A| = ^{l — cos6)d^(j) (defined in Eq. f l5^ ). A quantization 
condition for a monopole charge then arose naturally from the condition that a gauge field 
can be defined to cover the whole space (Eq. (!60|) ). In a case of spin damping monopole, an 
effective magnetic field in Eq. fl93|) cannot be written as a rotation of a local effective vector 
potential (because there is a monopole), but is represented by a nonlocal vector potential 
as Ss = V X A^""^^ + $, where (see Eq. ([8])) 

Af\r) = -P;Y^e,,J dr^N^ir - r'), (115) 

and $ represents a string singularity field. Based on this expression, the same argument 
as Dirac for a quantization condition [l| is valid, and thus the charge of a spin damping 
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monopole g is quantized as in Eq. (fT5|) . In order to understand fully the origin of spin 
damping monopole, it is important to study if the monopole is represented by patching two 
or more local gauge fields (see §VCp . This would be carried out by studying a clean (r = cx)) 
and strong sd coupling limit. 

A. Monopole interpretation of the inverse spin Hall effect 

The electromotive force discussed above is a result of the Maxwell's equation for the 
monopole, which is an exact equation required by the gauge invariance. If one detects 
the electromotive force given by Eq. f lll2p . it becomes a direct evidence of spin damping 
monopole. The electromotive force here acts on the spin of the electron, and thus can 



be detected electrically as has been demonstrated 15|. Remarkably, the observation of 



spintronics monopoles might have already been achieved. In fact, the electric voltage due to 



the magnetization precession has been observed in a system of ferromagnet on a Pt film [12 1 
and the direction of the voltage was in agreement with our prediction from the monopole. 
(The spin-orbit interaction expected in Pt is an intrinsic one due to periodic atoms, but 
the effect is expected to be the same as the one from the random potential. In addition, 
interface Rashba interaction might also be there. Therefore the monopole scenario would 



apply to the system in Ref. [12[.) In Ref. 12|. the mechanism for the voltage generation was 



argued to be the inverse spin Hall effect. According to the inverse spin Hall explanation, the 
magnetization precession generates spin current via the spin pumping effect, and the spin 
current, js, is converted into charge current by the spin-orbit interaction (the inverse of the 



spin Hall effect). This explanation assumes that the conversion mechanism of [44 1 



Ji = A5^e,,,4-, (116) 

where A is a constant representing the strength of the spin-orbit interaction and k is the index 
for the spin polarization of the spin current. This formula is, however, physically incorrect. 
In fact, the charge current is a conserved quantity, while the spin current is not because of 
the spin relaxation, and thus these two currents should not be simply proportional to each 
other. To put in other words, spin current does not have a unique definition because of its 
relaxation, and thus is not physical. A spin current explanation is therefore an approximate 
one which might be justified only at very short distance (less than the spin relaxation length). 
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Indeed, it has been theoretically demonstrated that that conversion formula is not satisfied 
in a case of slowly varying magnetization configuration 45| . In contrast to the spin current 
scenario, our monopole scenario is a result of Maxwell's equations, required by a symmetry 
of the electromagnetism. Since monopole is a conserved quantity satisfying Pm + V ■ jm = 0, 
monopole scenario is based on a relation between physical quantities, the electric field and 
monopole. Our scenario therefore explains the " inverse spin Hall effect" free from ambiguity, 
in contrast to the spin current one. 

For experimental confirmation of the spin damping monopole, of crucial importance is 
the separation of the monopole signal from the inverse spin Hall signal driven by the spin 
current. This is accomplished by applying an effective electric field perpendicular to the 
junction of Fig. [131 The monopole contribution then leads to the transverse electric current 
as a result of the Hall effect for the monopole, while the contribution of the spin current is 
not affected. 

Observation of the effective magnetic field emitted by monopoles, Eq. fl93|) . by use of 
muons or neutrons or electron holography would be also a direct evidence of monopole 
generation. 

B. Remarks on spin current 

Our analysis has also proved that the spin current does not modify the fundamental law 
of electromagnetism, namely Maxwell's equations. Historically, there have been arguments 
that the spin current would create the electromotive force by a modified Ampere's law based 
on the assumption that the spin current is equivalent to the fiow of two monopole charges 



461]. Quantum mechanics tells us that this idea is too naive, since the spin is a point object 
and thus the separation between the magnetic charges are strictly zero, even if one dares 
to interpret spin by two magnetic charges. The monopole current associated with the spin 
current is therefore absolutely zero. In fact, spin current is a quantity in an SU(2) spin 
space which couples to the electromagnetism only when properly projected, and since we 
know from the argument by Volovik 9|] that the projection in the adiabatic limit gives rise 
only to a hedgehog monopole current but not the spin current. Our study has proved that 
the correct projection of the magnetic systems with the spin-orbit interaction introduces a 
magnetic monopole in electromagnetism. 
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VIII. SUMMARY 



As we have discussed, magnetic monopoles are common objects in solids; they emerge in 
ferromagnetic metals from topological spin structures and from dynamic spin (magnetiza- 
tion) with damping. Although monopoles have been discussed so far based on a gauge theory, 
we have presented a novel method to access monopoles based on a transport calculation. 

In monopoles in ferromagnets, of particular interest are spin damping monopoles, gen- 
erated in a simple system of a ferromagnet and heavy metal. They open a novel path for 
connecting magnetism and electronics via an analog of Ampere's law. A novel concept of 
monopolotronics, control of monopoles, proposed in Ref. [lO|, is expected to be useful for 
realizing next generation spintronic devices. 



ACKNOWLEDGMENTS 

The authors are grateful to N. Kitazawa for discussion. This work was supported by 
a Grant-in-Aid for Scientific Research (B) (Grant No. 22340104) from Japan Society for 
the Promotion of Science and UK- Japanese Collaboration on Current-Driven Domain Wall 
Dynamics from JST. Two of the authors (A.T. and K.T.) are financially supported by the 
Japan Society for the Promotion of Science for Young Scientists. 

Appendix A: Electromagnetic fields in the relativistic representation 

In this section, we briefly summarize the electromagnetism theory in the relativistic repre- 
sentation using 4- vectors. We follow the convention used in Ref. [7| • Contravariant vectors, 
such as x^ = (t,x,y,z), are represented by the upper indices. Covariant vectors, repre- 
sented by lower indices, are defined by multiplying the metric tensor g^j^,^ as, for example, 
Xf, = g^.ux'', where 



g^L,. 



/l \ 
0-100 
0-10 

\0 0-1/ 



(Al) 
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Thus covariant vector of the 4- dimensional coordinate is thus x^ = {t,—x,—y,—z). The 
product of the covariant vector and contravariant vector is a Lorentz invariant scalar, e.g., 
Xfj^x'^ = t^ — r'^. (Note that we use relativistic notation only in this section. In other sections, 
upper and lower indices of the three-dimensional vectors means the same.) Differential 
operators are 



d.. 



d 



d d d d \ _ d 
dt^ dx' dy' dz J ' dx^, 



d d d d 
dt' dx' dy' dz 



(A2) 



We first describe the electromagnetism without monopole, namely, when there is. The 
electromagnetic field tensor is defined by a U(l) gauge field A^ as F^'^ = d^A'^ — d'^A^. 
Its components are (noting di = —d^) F^^ = —eijkB^ = Fij, where B'' = (V x A)^, and 



rOi 



{dfA + V0)* = — -E* = — -Foi- III the matrix representation. 



Ff"" = d^'A" - d^A^" 



( 



-Er,. £/„ E, ' 







By B^ 



-Bz By 

. -B^ 
\Ez -By B, / 

The field strength tensor satisfies by definition the following identity 



(A3) 



d^F^" = 0, 



(A4) 



where (6*^123 = j^ = —60123 is the four- dimensional antisymmetric tensor) 
1 



ptii/ ^ Z^t^'^P^ p 



pa- 



(A5) 



The 0-component of the identity (IA4p is 



5.F* = — e'^^ViF^^ = V ■ S = 0, 

2 



and i-component reads 



9oF°^ + d,F^' = - {doe''''F,k - 2dje'^'Fok) 

d 
= -—B-V X E = 0. 

dt 



(A6) 



(A7) 



Therefore, the condition of no monopole, V ■ -B = 0, and the Faraday's law are trivial result 
of the U(l) gauge symmetry (definition of F'^'^ = d^A" — d" Af^). 
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Appendix B: Spin damping 

Let us look into a role of spin damping represented by a vector N = S x S. A spin 
dynamics is induced by magnetic field. An equation of motion for a spin is thus generally 
given by Landau-Lifshitz equation 

S = 7Stot X 5, (Bl) 

where 7 is gyromagnetic ratio and B^ot represents the total magnetic field acting on S. 
Btot thus includes an external magnetic field, an internal field due to exchange interaction 
with other localized spins, magnetic anisotropy field, and also the effect of coupling to other 
degrees of freedom such as conduction electrons and phonons. 

We know that damping or dissipation arises in general from a coupling to an environ- 
ment, namely, other degrees of freedom. In a case of spins in metallic magnets, most strong 
source of dissipation is conduction electron. Dissipation in this case is caused by spin-orbit 
interaction, which converts spin angular momentum into orbital angular momentum. There- 
fore, damping effect is calculated by evaluating an effective magnetic field from conduction 
electrons, resulting in |47| 

7Sei = -aS, (B2) 

where a is a constant proportional to the spin-orbit interaction. The damping torque is 
thus represented by aS x S, called the Gilbert damping term. The equation of motion thus 
reads 

S = -fBx S + aSx S, (B3) 

where B = Btot — -Bcnv is the field neglecting the effect of the environment, -Bcnv Equation 
( ]B3l) is called Landau-Lifshitz-Gilbert equation. 

The effect of Gilbert damping torque is understood by looking into Fig. [121 In fact, a 
vector S X S = N tends to point the spin perpendicular to its precession direction, i.e., to 
the equihbrium direction along the field B, and hence aN represents the spin dissipated by 
the environment. 
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